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Abstract 

A huge swath of parameter space in the context of the Minimal Supersymmetric Standard 
Model (MSSM) has been ruled at after run I of the LHC. Various exclusion contours in the 
rrig — m^o plane were derived by the experimental collaborations, all based on three-body 
gluino decay topologies. These limits are however extremely model dependent and do not 
always reflect the level of exclusion. If the gluino-neutralino spectrum is compressed, then 
the current mass limits can be drastically reduced. In such situations, the radiative decay of 
the gluino g —)■ gxi can be dominant and used as a sensitive probe of small mass splittings. 
We examine the sensitivity of constraints of some Run I experimental searches on this decay 
after recasting them within the MadAnalysisS framework. The recasted searches are now 
part of the MadAnalysisS Public Analysis Database. We also design a dedicated search 
strategy and investigate its prospects to uncover this decay mode of the gluino at run II 
of the LHC. We emphasize that a multijet search strategy may be more sensitive than a 
monojet one, even in the case of very small mass differences. 


PACS numbers: ll.SOpb, 14.80Ly, 14.80Nb 
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1 Introduction and motivation 

During run I of the LHC, in addition to the quest for the Higgs boson, an extensive search 
programme for New Physics (NP) phenomena has been carried out. However, no sign of 
physics beyond the Standard Model (BSM) has been unravelled so far [1-4]. As a result the 
parameter space for weak scale NP, including models of weak scale Supersymmetry (SUSY), 
has shrunk considerably. 

The simplest incarnations of weak scale SUSY, the MSSM and its variants, depending on the 
SUSY breaking scheme, are facing a wealth of experimental data [1-4], including the Higgs 
boson mass measurement mu — 125 GeV [5]. Since no NP signal was revealed, the only 
appropriate physical interpretation of these searches was to set limits on cross sections for 
the production of SUSY particles and then derive bounds on their masses. Interpreting the 
NP experimental searches is a non-trivial task and almost impossible to perform in a model- 
independent way. The early searches for SUSY phenomena have focussed on the constrained 
MSSM (CMSSM)/minimal supergravity (mSUGRA) and probed mainly the gluino {g) and 
the squarks (g) of the hrst two generations [6,7]. The current lower limits on the gluino and 
the hrst two generations squark masses in the framework of the GMSSM stand at rrig > 1.7 
TeV for almost degenerate gluino and squarks and rriq >1.4 TeV for very high squark masses. 
However, the universal boundary conditions imposed at the Grand Unihed Theory (GUT) 
scale in such models impose particular relations between particle masses, decay branching 
ratios, etc... which severely restrict the way a SUSY signal could show up. To relax some of 
these assumptions, the Simplihed Model Spectra (SMS) [8-13] approach has been adopted 
systematically by the ATLAS and GMS collaborations to interpret the results. Instead of 
considering the full SUSY spectrum, a SMS search only targets one or two relevant decay 
modes (and the particles involved in) with 100% branching ratio (HR). In turn, a handful of 
parameters are to be considered for interpreting the analysis: the production cross section 
of the mother particle, its mass, and the ones of the daughters. Nevertheless, such SMS 
approaches may not cover all patterns in which a NP signal could be unveiled, especially 
if one consider the plethora of existing supersymmetric models. More precisely, even in R- 
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parity conserving scenarios, if the lightest SUSY particle (LSP) is massive and degenerate 
with the squarks and/or gluinos, the so-called “compressed” SUSY scenarios, these limits 
can be seriously weakened [14-20]. 

As far as the official experimental analyses are concerned, published gluinos mass limits in 
the SMS approach were obtained mainly from its three-body decays g —)■ ggy?, where the 
hnal states quark can be of any flavour [6,21-25] or ^ qq'xf [6,21]. Provided the gluino- 
neutralino spectrum is compressed, the limits derived from three-body decays involving 
heavy-flavours do not apply and only the ones with light-flavours in the hnal states can 
give a handle on this situation. However, as limits on the hrst and second generation of 
squarks are quite strong, the g —)■ qqxi modes with light quarks may be heavily suppressed. 
Furthermore, as already stressed, a large part of the parameter space feature a compressed 
spectra to evade the stringent bounds from the LHC, and few SMS interpretations fully 
encompass this situation for the gluino. In compressed scenarios, initial state radiation 
(ISR) in the form of jets which recoil against the missing momentum of the LSP may be 
used to discover or place bounds on the model. This method has been used to set mass 
limits on squarks and gluinos [6,14-20,26-30] and third generation squarks [14,29,31-35]. 
In the case where the three body decays are kinematically forbidden or suppressed, there 
exists another decay pattern available for the gluinos to decay into, which has received very 
little attention so far, but could £11 the small mass difference gap: that of the two-body 
radiative decay g gXi- 

In a CMSSM-inspired scenario, since the mass difference AM = rrig — is fixed from con¬ 
ditions imposed at the GUT scale, when rrig increases, so does AM. In turn, the decay width 
of the gluino will be dominated by the three-body decays. Unlike GUT constrained SUSY- 
scenarios, general explorations of the MSSM may exhibit significant portions of parameter 
space where the radiative decay g —)■ gxi can be competitive with three-body decays [36-45]. 
Within these specific regions, significant branching fractions of the gluino radiative decay are 
characterised by a spectrum with heavy squarks. Thus, the direct two-body decays g ^ qq 
are forbidden and g qqx^ are suppressed by the heavy scalar quark mass, leaving the 
loop-induced decay as the only available mode to decay into. Therefore, for compressed 
gluino-neutralino scenarios, the gluino loop decay can be used as a sensitive probe and 
improve the current limits. 

A spectrum with decoupled scalars and lighter gauginos and higgsinos is not unexpected in 
SUSY models like (mini) Split-SUSY/PeV scale SUSY [46-52], Spread SUSY scenarios [53] 
or pure gravity mediation [52,54-56]. Moreover, a Higgs mass around 125 GeV implies 
a non-negligible branching fraction of the gluino loop decay when tau/d ~ C>{1) [43]. It 
is worth mentioning that if the scalar superpartners are heavy enough, then the gluino 
lifetime Tg can be sufficiently large to lead to displaced vertices or even hadronise into R- 
hadrons [49,50,57,58]. Nevertheless, if m < 0(10^) TeV, Tg is too small such that the gluino 
decays promptly in the detector [43]. 

To obtain a compressed gluino-neutralino spectrum, given a SUSY-breaking scheme, one 
needs non-universalities in the gaugino sector at the GUT scale to obtain a small AM. In the 
context of non-universal supergravity (NUSUGRA) models, the gluino can be in a significant 
part of the parameter space the next-to-lightest SUSY particle (GNLSP) and the neutralino 
the LSP [59-66]. A non-exhaustive list of other SUSY models accommodating the GNLSP 
scenario include left-right symmetric model with gravity mediated SUSY breaking [67-70], 
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pure gravity mediation/mini-split SUSY [52,56] and general gauge mediation [71]. 

Besides being a discovery mode for the gluino, its loop-decay has some other merits from 
the phenomenological perspective. Indeed, since squarks contribute to the radiative decay 
through loop diagrams, the two-body decays g —?■ gXi carry important information about the 
scalar mass scale fh [43], which would be otherwise inaccessible in scenarios with decoupled 
scalars. Moreover, in GNLSP models, the LSP can account for the relic density of dark matter 
through its coannihilation with the gluino. The efficiency of such processes crucially depends 
on the gluino-LSP mass difference AM. Thus, measuring or constraining this splitting, only 
accessible at collider experiments, is essential to estimate the phenomenological viability of 
these models. 

We should also note in passing, that in models with spontaneously broken global supersym¬ 
metry the radiative decay g gG, where G is the gravitino, is also possible. It can be 
dominant when the original scale of supersymmetry breaking y/F is of the same order of the 
squark mass scale rh [42,48,72-79], whatever the rrig — ttIq mass difference. Depending on 
the value of F the gluino can decay promptly, in flight inside the detector, or be long-lived. 
For very long-lived gluinos, ATLAS and CMS have excluded charged R-hadrons with a mass 
less than 1.3 TeV [80,81], if they escape the detector, or less than 850 GeV [57] if stopped 
in the detector. If metastable, the gluino is excluded up to rrig = 850 GeV for decays to 
qqXi/9Xii for ^ lifetime of 1 ns and m^o = 100 GeV [82]. 

A work probing scenarios where the LSP and the gluino are almost degenerate was recently 
performed in [29]. However, the limits derived there solely considered the three-body decays 
of the gluino into light-flavour quarks g qqXi: but with light sbottoms and light stops in 
the spectrum. As we will review in the forthcoming sections, the radiative decay g —?■ gxj is 
favoured when light stops are present and should dominate instead. Therefore, investigating 
rather the gluino loop-induced decay seems better motivated from a phenomenological point 
of view when all scalars are heavy and AM small. The ATLAS collaboration has also recently 
provided, using only a hadronic search that we reimplemented for this work, some results 
on the radiative gluino decay in [6]. Yet, results have only been presented for a hxed gluino 
mass rrig or hxed LSP mass rn^o, but not in the full rrig — m^o plane. 

In this work, we aim at quantifying to which extent the radiative gluino decay can constrain, 
using Run I monojet and multijets analyses, the region where the mass splitting AM = 
rrig — m^o is small, assuming a SMS interpretation where BR(^ fi'X?) = 100%. We argue 
in the next section, after investigating the parameter space where the loop decay is important, 
why a somewhat lighter third generation squark spectrum is desirable to enhance this decay 
from the phenomenological MSSM (pMSSM) point of view. We use the MadAnalysis5 [83,84] 
framework to derive exclusion contours in the rrig — m^o plane, deduced from a Simplihed 
Model approach, after recasting ATLAS and GMS monojet and multijets analyses. We 
provide in our work a full quantitative exploration of this plane and the quoted ATLAS 95% 
GL upper bound [6] for hxed rrig and m^o rather serves to us as a cross check on our results, 
which give us quite conhdence in our procedure with which we derived our own mass limits. 
These recasted analyses are now available on the MadAnalysis5 Public Analysis Database 
(PAD) [85]. The last section is devoted to the determination of a search strategy at 13 TeV 
to investigate the prospects for discovering the gluino through its radiative decay at run II 
of the LHG and then we conclude. 
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2 Phenomenology of the radiative gluino decay 


In this section we review the gluino decays in a heavy scalar scenario. When the squarks are 
heavier than the gluino, the latter can undergo three body decays into a pair of quarks and 
a neutralino or chargino through an off-shell intermediate squark, or decay radiatively into 
a gluon and a neutralino, 

g^q + q + X^i ( 1 ) 

q + q' + xf (2) 

^ + X? (3) 

where [i = 1... 4) and xf U = 1)2) are neutralinos and charginos respectively. Let 
us first review the radiative two body decays of the gluino, which occm^via quark-squark 
loop diagrams (Fig. 1). The neutralino is an admixture of wino- (IPs), bino- {B) and 



Figure 1: Diagrams relevant for the radiative gluino decay. Since the gluino is a Majorana 
fermion each diagram is accompanied with a similar one but with the internal arrows reversed. 


higgsino-like {H) neutral spinors, which are respectively the superpartners of the neutral 
gauge eigenstates B, IF 3 and of the two Higgs doublet of the MSSM Hu d, 

+ W2IP3 + Ni,Hd + N^aHu ( 4 ) 


If the squarks are very heavy, from an effective Lagrangian point of view the decay width into 
a wino is strongly suppressed since it is induced by a dimension 7 operator [42,43], whereas for 
a bino or a higgsino it can be induced already by a dimension 5 chromo-magnetic interaction, 

CeE. = (5) 

m ^ 

where fh is an effective squark mass scale. The partial widths of these two body-decays are 
given by [36,42,43,69], in the heavy squark limit. 
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where a and Og are the electromagnetic and strong couplings constants respectively. We 
used the short-hand notations cw = cos 9w,sw = sin 6 *w,S /3 = siu/S. The angle 9w is 
the weak mixing angle and f3 is defined through the ratio tau/S = Vu/vd where Vu,d are 
the vacuum expectation values of the two SU{2)l Higgs doublets of the MSSM. Both 
angles lie between 0 and 7 r/ 2 . Yq, appearing in Eq. 6 , is the quarks hypercharge where the 
index q runs through all flavours. In Eq. 7, since the decay width is proportional to the 
top mass, only the stops-tops loops were kept. From Eq. 6 we can see that r(^ —?■ gB) 
scales like whereas r(^ —?■ gH) benehts from a logarithmic enhancement. The origin of 
this logarithmic enhancement can be understood from an effective Lagrangian point of view 
as follows. In the effective theory where the squarks have been integrated out, this decay 
is induced by a top-top-gluino-higgsino operator in which the two top quarks form a loop 
which can emit a gluon to form the chromo-magnetic gluino-gluon-higgsino interaction of 
Eq.5. Such a diagram is divergent in the effective theory thereby generating a logarithmic 
enhancement [36,42], which is cut-off by the squark mass (the scale of the effective theory 
breakdown). Therefore this decay is an interesting probe of the scalar mass scale in case 
they are kinematically inaccessible from colliders [43]. In the very heavy squark limit the 
leading logarithmic corrections should be resummed to obtain reliable predictions using 
renormalisation-group techniques as has been done in [42]. 

Let us now turn on to the three body decays of the gluino g —?■ qqx^ which is mediated 
by squark exchange q (Fig. 2), with the assumption that the squark sector respects flavour 
symmetries. In the massless quark limit, only the bino-like contributes and we can get an 



Q 



Figure 2: Three-body decay of the gluino where the (s)quark flavour runs through q = 
M, d, c, s, f, b. 


analytical expression of the gluino partial decay width [39,42,43,69], 


r(^ ^ qiq'kB) 


aasYq^ ml 

9671 cos^ 9w 



( 8 ) 


where f{x) = 1 — 8 a; — 12a;^ In x -|- 8 x^ — and g{x) = 1 -|- 9x -|- 6 x In x — 9x^ -|- 6 x^ In x — x^. 
Like the two-body decay into the bino-like neutralino in Eq. ( 6 ), the three-body decay width 
scales as in the heavy squark limit. If the first two squarks generations are decoupled 
from the third one, which is quite easily achieved in SUSY since the large top Yukawa usually 
drives the running of third generation squark masses to lighter value than the two others, the 
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quarks produced in the three body decays will be composed of tt, bb pairs. For generic quark 
masses the integration over the three-body phase space cannot be performed analytically 
and one must resort to a numerical integration. However, keeping only the quark mass 
dependence in the couplings and in the limits m^o -C rrig and rrig <C the decay 

rate g —)■ qqXi [37] is given by 


r(^ ^ qqX°) ^ 
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The sum is over all quark species allowed by kinematics. The a, b, c coefficients are given by 

gmuNi4^ grudNis 


^7/ 


2Mw sin (3 ’ 2Mw cos (3 ’ 

bqL = 9 fan OwNa + , bg^ = gQq tan OwNn 


( 10 ) 

( 11 ) 


where u = u,c,t and d = d,s,b. The electric charge Qq is defined as Qq = -|- Yq/2. 

It is then natural to ask oneself in which cases the two-body decay can dominate over the 
three-body one. To investigate this we define the ratio R 2/3 as 


D _ ^{9 ^9X°) 

r(^ ^ qqX^) 


( 12 ) 


Obviously, the three-body decay is suppressed by the phase space when the mass difference 
between the gluino and the neutralino is too small. This is particularly true for g —)■ 
when the tt threshold is closed, i.e rrig — m^g < 2mt. In this situation the decay rates for 
g —)■ bbx^ may still compete with r(^ ^ 9 X^)- In the higgisno-like case, the bb and gx^ final 
states are dominated by the heaviest quark: the former is proportional to ml and the latter 
to mj. In the limit where the left and right squark masses are equal for each respective 
flavour q, i.e mq^ = = fhq, using Eqs. (7) and (9) and defining xt = ml/ml, one finds 

for R 2/3 [37] 
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(13) 


In the bino-Iike case, we would obtain a similar formula using Eqs. (6) and (9), but R 2/3 
only benehts from the parametric enhancement rhb/mt. Hence, there are two conditions for 
which the two-body decay will dominate over the three-body one: 


• If the tt threshold is closed one needs a somewhat hierarchy between the sbottoms and 
the stops, 

• Decouple all the squarks at very high masses, in this case R 2/3 oc mf/m?[l -|- Inx^]^. 

However, as stressed in [42], for such very heavy squarks, the logarithmic corrections have 
to be resummed and in the case of all squarks masses being equal and large, the resum¬ 
mation tends to decrease R 2/3 compared to Eq. (13) when the common scalar mass scale 
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fh > 10® GeV. On the other hand, in the bino-like case, i? 2/3 does not benefit from such 
parametric enhancements. The lack of logarithmic enhancement for the two-body decay {i.e 
the higgsino-like neutralino is not allowed in gluino decays) implies that the three-body will 
generally dominate instead, if available. 

It is worth adding that in the pure higgsino case, it is rather another three-body decay which 
rival with the gluino loop decay, namely g —)■ tb/tbxf- This is due to the fact that the neutral 
and charged higgsinos are almost degenerated in mass, and this decay will dominate only 
once the th/ib threshold is crossed, i.e rrig — m^o ~ rrig — Tn^± > mt + mi,. These three-body 
decays do not suffer from a hierarchical third generation squark spectrum since they are 
mediated by stops. In the pure bino case this decay will be effective at larger AM since the 
lightest chargino Xi is much heavier than the bino-like LSP. 

In order to illustrate to which extent the radiative gluino decay can dominate over the 
other decay modes, without assuming any specific SUSY-breaking pattern, we performed a 
numerical scan within the pMSSM parameter space using SUSY-HIT [ 86 ]. In this numerical 
investigation we set all soft masses (sleptons and squarks) to 5 TeV and trilinear couplings 
Ai to zero except for the third generation squarks for which we chose the following range: 

• = 1 TeV, Mq^ = 2 TeV, Ab = 0 

• 1 TeV < < 2 TeV 

• —2 TeV < At < 2 TeV 

The parameters related to the gaugino and Higgs sector were set to 

• 400 GeV < Ml, fi < 800 GeV 

• M 2 = 2 TeV, M 3 = 600 GeV, tan /3 = 10 

The scalar mass scale chosen is sufficiently low such that we ensure that the gluino would 
decay in the detector. We picked up only points giving the SM-like Higgs mass in the range 
122 GeV < mh < 128 GeV and the resulting gluino mass is driven around ~ 810 GeV 
due to radiative corrections. We ensured that > m^o such that the lightest neutralino 
is always the LSP. We divided the surviving points of the scan into three different regions, 
depending on the nature of the lightest neutralino, whether it is mostly higgsino-like, bino-like 
or a mixed bino-higgsino. The higgsino-like (bino-like) neutralinos correspond to a higgsino 
(bino) fraction greater than 90%. The remaining points are an admixture of bino-higgsino 
like neutralinos which we dubbed as “mixed”. 

We display in Fig. 3 the branching ratio of the two body decay of the gluino into a gluon plus 
the lightest supersymmetric particle (LSP) in terms of the mass difference AM = —m^o. 

Given the nature of the spectrum, only the radiative decays g —)■ gXi and the three-body 
modes into heavy quarks {g —)■ bb/tixi g —)■ tb/ibxi) are allowed. The colour shading 
correspond to the level of hierarchy between the lightest stop ti and lightest sbottom bi 
given by the ratio It is worth noting that in the higgsino-like case we also included 

the two-body decay of the gluino to the next-to-lightest neutralino g —)■ gx 2 since the mass 
difference between the Xi and y® is at most ~ 15 GeV. For such a small mass difference the 
next-to-lightest neutralino y® would decay into the LSP and soft leptons or quarks originating 
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Figure 3: Branching ratio of g ^ gx^ in terms of AM = rUg — m^o for a higgsino-like 
(top,left), bino-like (top,right) and mixed bino-higgsino (bottom) neutralino. In the higgsino- 
like case the branching ratio of the radiative deeay g —?■ gx 2 'Is also included. The eolour 
shading indicates the ratio 


from the off-shell Z boson decay in yf) —)■ Z*^. Hence the channel g —)■ gx^ would leave 
the same experimental signature of jets -|- MET in the detector as ^ gxfl. Glancing at 
Fig. 3 it is clear that largest branching fractions for the radiative decay are obtained when 
the mass difference AM is small. As one could expect from the approximate formulas in 
Eqs. (6),(7), the branching ratio of the two body decays g —)■ gHi 2 (top, left panel of Fig. 3) 
gives the largest branching fractions. This is due to the fact that it depends on the fourth 
power of the top mass, unlike the bino-like and mixed cases, see Eqs. (6),(7). In addition, 
the dependence on ratio of the loop-induced decay is much more pronounced in the 

bino-like and mixed case than in the higgsino case. It is particularly true for the bino case 
(top, right panel) where largest branching fraction are only obtained when the hierarchy 
between the sbottoms and stops is signihcant. Moreover, the two-body decay into higgsino- 
like neutralinos dominate even for quite large AM: the radiative decay starts to fall off at 
AM ~ 200 GeV, while for the bino-like and mixed cases they already decline at AM ~ 100 
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GeV. The drops of the loop-induced decays can be understood from looking at Fig. 4. In 
this hgure are displayed the branching fractions of the most important three-body decays 
for each nature of the LSP. The left panel features the branching ratio BR(^ bbxi) for 
bino-like or mixed neutralinos. The higgsino-like LSP case is not shown since it is negligible 
and buried in the two others. The right panel presents the branching fraction of the gluino 
decay into th/ib final states for all three LSP natures considered. In the higgsino case (blue 
diamonds in the right panel of Fig. 4), a comparison with the top, left panel of Fig. 3 makes 
it clear that as soon as the decay into th/ib is kinematically allowed, the branching fraction 
BR(^ —)■ gxi 2 ) decreases fast (for AM > 200 GeV) and the three-body mode g —)■ tb/ibxf 
starts to take over. However, thanks to the logarithmic enhancement, the loop-decay of the 
gluino to neutral higgsinos still persists at large AM and is of order 10 — 20%. The rationale 



’ Higgsino LSP 
■ Mixed LSP 


0.8 — ■ Bino LSP 



AM [GeV] 


Figure 4: Branching ratios for the most important three-body decay modes depending on 
the composition of the LSP. The left panel corresponds to the channel g bh/d, where 
red squares correspond to the case where the LSP is bino-like and green stars when it is 
a mixed bino-higgsino. The pure higgsino case is not shown as it is negligible. The right 
panel displays the g —?■ tb/tbxf for each composition, blue diamonds points corresponds 
to situations where the LSP is higgsino-like, orange points and when it is mixed magenta 
squares. 

for this decay mode to take over in the higgsino case is that, as soon as the threshold is open, 
it does not suffer from large suppression factors. Indeed, this decay is also mediated by stop 
exchange, as the loop decay is, and not suppressed by small quark masses. In the pure bino 
case, BR(^ —)■ tb/ibxf) is suppressed until larger AM. This is due to the parameter range 
we chose for the scan. Indeed, for this decay mode to be open one needs a significant mass 
splitting between the gluino and the lightest chargino, typically — m-± f mt -\- mf,. In 
our scan m-± is driven by the value of the parameter /i, which must still be quite large to 
keep the LSP bino-like. Therefore one needs to reach a significantly larger mass splitting 
AM to fulhl these requirements, compared to the pure higgsino case where m^o ~ m-± ~ u. 
Instead, the three-body decay mode g —> bhxi quickly dominates over g —)■ gx^, already from 
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AM ~ 100 GeV. This occurs even when the ratio is large. Again, the underlying 

argument is that the gluino loop decay to a bino-like neutralino does not proht from any 
enhancement, apart from phase space considerations. For a mixed neutralino, when AM 
increases, the suppression of the radiative decay is counter-balanced by both g —)■ hhxi and 
g tb/tbxf. 

Having identihed in which region of parameter space of the pMSSM the radiative decay can 
dominate, we now turn on to which extent it can be constrained by existing ATLAS and 
CMS 8 TeV analyses. 

3 Constraining the gluino loop decay from existing 8 
TeV LHC analyses 

It is expected that the typical signature of the radiative gluino decay will lead to a prepon¬ 
derance of jets. If favoured, when the mass difference AM = rrig — m^o is small, the hnal 
state jets are typically as soft as those originating from the parton shower and in turn the 
gluino decay gets buried into the huge QCD background. To increase the sensitivity one can 
resort to, similarly to what is done concerning Dark Matter or compressed spectra searches, 
trigger on a hard initial-state-radiation (ISR) jet, the so-called monojet searches. 

Since the radiative gluino decay has not been fully investigated at the LHC, we wish to hll 
this gap, using in particular monojet analyses to be sensitive to it. It is worth to note that 
the multijet analysis performed in [6] has indeed probed a small mass AM. Nevertheless, 
the interpretation has been done in a different context from ours, where the hrst and second 
generation squarks are not too heavy and therefore the three-body decay into light-flavour 
quarks g —)■ qqxi is not suppressed. Thus, one can probe small mass splitting between the 
gluino and the LSP in this case. In [6], the direct decay of the gluino to a gluon and a 
neutralino was also used to constrain rrig and m^o, but based only on a fully hadronic search. 
Only numbers for particular values of rrig and m^o were presented. For a massless neutralino, 
gluino masses below 1250 GeV can be excluded, and a lower mass limit of m^o = 550 GeV 
corresponding to rrig = 850 GeV is quoted, but no full coverage of the rrig — m^o mass 
plane has been presented. In [29], an exclusion contour in the rrig — m^o mass plane has 
been derived using monojet searches, using the process g —)■ qqXi assuming 100% branching 
fraction, among other scenarios investigated. However, as we saw in the previous section, the 
dominance of this channel strongly depends on the squark spectrum, which is not detailed 
in [29], except that the stops and sbottoms are assumed to be light. We saw in the previous 
section that the radiative decay is favoured when light stops are present and therefore it is 
not clear if the assumption that g —)■ qqxi dominates holds. 

In the context where the hrst and second generation of squarks are degenerate, bounds 
on their masses are pushed beyond the TeV scale. In the scenario we are focussing on, 
namely with decoupled hrst and second generations with somewhat lighter third generation, 
this ehort is complementary to the official experimental analyses and worth to investigate. 
Achieving such a task then requires to recast the already existing ATLAS and CMS analyses 
since we do not have access to the raw data of both collaborations. The experimental 
signature for the gluino radiative decay is rather simple: it only consists of jets (possibly 
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very soft) and missing transverse energy (MET). A characteristic feature of scenarios where 
the two-body radiative gluino decay dominate is typically the ones which are difhcult to probe 
in a collider environment, namely ones where the mass difference between the mother and 
daughter particles is small (he compressed spectra). These scenarios are difficult to explore 
since they lead to hnal states with less energetic states or leptons, thus reducing the detection 
efficiency and signal acceptance. This may be part of the reasons why this decay pattern has 
not been fully considered as worth to investigate hrst in the past, as a priori its prospects are 
not very enthusiasming. In addition, in typical mSUGRA/CMSSM scenarios, this branching 
ratio is quite small since the mass splitting between the gluino of the LSP is generically 
quite large. However, the advent of monojet searches to precisely increase the sensitivity to 
compressed spectra has been a game changer in this respect. As noted in [43], when the ratio 
of the Higgs doublet vev’s tan f3 ~ 0{1), a 125 GeV Higgs implies a non-negligible branching 
fraction of the radiative modes for the gluino decay. Moreover it is well-motivated in SUSY 
scenarios where the gluino is the next-to-lightest supersymmetric particle (NLSP), the so- 
called GNLSP models. Examining this process at the LHG can also put constraints on the 
viability of these models in solving the Dark Matter problem. Indeed, in such models, the 
right amount of relic density is obtained through important coannihilation processes between 
the LSP and the gluino and their efficiency is mainly driven by AM see for example, [87] 
and references therein. Hence, a full quantitative investigation of the radiative gluino at the 
LHG seems timely. 

The experimental signature of the radiative decay can be categorised in 1, 2 and 3 jets bins 
plus MET. Unfortunately there exists no published official analysis covering all such cate¬ 
gories. In this purpose, and to cover all jet bins, we reimplemented within the MadAnalysisS 
[83,84] framework the following three analyses: 

• ATLAS-SUSY-2013-21 [88]: this analysis targets the stop decay ti c + Xi using a 
monojet analysis. It has also been used to probe the small squark-neutralino mass 
difference using q ^ q + Xi^ see [6], 

• ATLAS-SUSY-2013-02 [89]: search for squarks and gluinos in hnal states using high-p-r 
jets with no leptons, 2-6 jets and large MET, 

• CMS-SUS-13-012 [90]: search for new physics in multijet events with large MET di¬ 
vided into three jet multiplicity categories (3-5, 6-7, and 8 jets), 

in order to work out the 8 TeV LHG constraints using 20.3 fb“^ of collected data. The 
recasted GMS-SUS-13-012 SUSY search was already available in the MadAnalysisS Public 
Analysis Database [85] and can be found at [91]. The two remaining searches were im¬ 
plemented and validated for the present work and are now available at [92,93]. For the 
sake of clarity, we present in the Appendix A a comparison between the official cut hows 
and our reimplementations of the ATLAS monojet [88] and multijet [89] analyses within 
MadAnalysisS (MAS). From Tables 4,7 and 8 the selection cuts of each respective analysis 
can be read oh. 

For the monojet search, see Table 4, one can see that the agreement is not very good at 
the level of the preselection cuts. The discrepancy originates already from the E™'®® > 100 
GeV cut applied at the reconstructed level. The MET trigger efficiency varies with U™'®® 
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and in particular is not 100% effective at = 100 GeV. To reproduce the MET trigger 
efficiency we parametrised the efficiency turn-on curve presented in [94] coming from the 
ATLAS simulation of the process pp ^ ZH ^ uubb, as advised after communication with 
the ATLAS SUSY conveners. We already observe at this level a discrepancy which ranges 
from 16% to 32% for the three benchmark points between our implementation and the 
ATLAS result. Moreover we cannot reproduce the Trigger, Event Cleaning, and Bad Jet 
veto efficiencies, as we do not have access to this information in official documentation. 
However, once we enter the signal regions, the relative efficiencies of our cuts are comparable 
with the official ones and at the level of hnal number of events we observe an agreement 
between 15-25%, which is in the ballpark of expected accuracies from fast-simulation. More 
details about the validation of the recasting procedure can be found at [95]. 

For the ATLAS multijet analysis [89] a cut flow comparison is also presented in Tables 7 
and 8. One can see that the agreement between the recasted analysis and the official is 
fairly good. The CMS multijet implementation and validation was already available on the 
PAD [85] before the initiation of this work and all details were presented there. 

As official results are not provided in terms of the gluino branching ratios and the radiative 
mode has not been fully investigated so far, we derived limits on this decay pattern within a 
Simplihed Model approach. In the case we focussed only the gluino and LSP are relevant in 
the mass spectrum and we assume BR(^ —)■ gxi) = 100%. The signal generation was done 
using MADGRAPH5 [96,97] using the CTEQQ6L1 parton distribution functions [98], as 

+ + (14) 

where X = 0,1 is the number of generated jets in the signal sample. The MADGRAPH samples 
are produced using the AUET2B tune [99,100]^. The MLM [102] matching scheme is used, 
and a MADGRAPH kx measure cut-off and a PYTHIA [103] jet measure cut-off both set to 0.25 
times the mass scale of the SUSY particles produced in the hard process. Hadronisation is 
performed using PYTHIA-6.426 [103] as incorporated in the MADGRAPH-5 package and detec¬ 
tor effects are simulated with the ’MA5tune’ of DELPHES-3 [104] as described in [84]. The 
number of events are then normalised to the correct luminosity after including cross sections 
at the next-to-leading order and next-to-leading logarithm accuracy [105], as tabulated by 
the LHC SUSY Cross Section Working Group [106] 

For the statistical interpretation, we make use of the module exclusion_CLs .py provided 
in the MadAnalysis5 recasting tools. Given the number of signal, expected and observed 
background events, together with the background uncertainty (both directly taken from the 
experimental publications), exclusion_CLs .py determines the most sensitive signal region 
(SR), the exclusion conhdence level using the CLs prescription from the most sensitive SR, 

^It is worth to note that care should be taken when using recent versions of MadGraphS together with 
the AUET2B tune. Indeed this tune was designed from a Leading-Order (LO) Matrix Element -I- Parton 
Shower (MEPS) merging. Since MadGraph5.2, the matching MEPS is done at next-to-leading order and in 
the same fashion as in MCNLO [101]. Therefore if one uses this tune together with a version 5.2 or higher of 
MadGraph, the interpretation of the results should be taken with great care. Indeed, we observed signihcant 
discrepancies in reproducing the official exclusion curves between using a recent version of MadGraph or the 
one actually adopted in the official analyses in [88,89]. That is why, for each analysis, we stick as much as 
possible to the Monte-Carlo conhgurations described in the experimental analyses. Thus, each signal sample 
has been generated with a different MadGraph version for each analysis we recasted. 
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and the nominal cross section cxgs that is excluded at 95% CL. The simplihed model parameter 
space is simply scanned by varying rrig and m^o until a 95% CL excluded number of BSM 
signal events is reached. The derived limits in the {mg,m^o) plane are presented in Fig. 5. 
In the left panel of Fig. 5 the red lines are the 95% CL exclusion bounds using the recasted 



Figure 5: Left panel: 95% CL exclusion contours for the radiative gluino decay simplified 
topology. The solid red line corresponds to the mass limits obtained from the MAS recasted 
ATLAS monojet analysis [93], the red broken line from the MAS recasted ATLAS multijet 
search [92] and the dashed-dotted line from the MAS recasted CMS multijet analysis ]91]. For 
comparison, the official 95% CL exclusion lines for the g tix\ SMS from ATLAS (green 
dashed line) ]6] and CMS ]23](blue dashed line) as well as for the g QQXi SMS from 
ATLAS ]6](purple dashed line) are also shown. Right panel: 95% CL exclusion contours in 
the plane mg versus AM where AM = mg — m^o. Only the exclusion contours from the MAS 
recasted monojet and multijet analyses are given since they provide the best sensitivity in the 
low AM region. 

analysis^. The results of the three recasted analyses (solid line : ATLAS monojet [93], 
broken line: ATLAS multijet [92] and dashed-dotted: CMS multijet [91]) are compared 
against three official published results: in green and blue are displayed the 95% CL limits 
obtained from the simplihed topology g —)■ tix] from ATLAS [6] for the former and from 
CMS for the latter [23]. One immediately observes that the official SMS interpretations are 
limited when the tf threshold of the gluino three-body decay is closed. The CMS results 
pushes the limit a little bit above the tf threshold limit using off-shell tops, as compared 
to ATLAS. The purple broken line corresponds to the ATLAS limits in the {mg,m^o) mass 
plane for the g —?■ qqx] simplihed topology using the multijet search [6]. In the corresponding 
model investigated all squarks are decoupled and since the hnal states quarks are massless in 
this simplihed approach, a small AM can be probed. In our scenario only the light-havour 

^Note that one could have used, for the monojet analysis, the acceptance times efficiency tables as provided 
on the Twiki page of [88]. We chose to recast anyway this analysis to enrich the PAD [85]. 
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squarks are completely decoupled from the low-energy spectrum and therefore our results 
are complementary to these. 

A word of caution is in order before discussing and interpreting our results. Indeed, we do 
not provide an estimation of the theoretical uncertainties on both the signal cross section 
and acceptance. This certainly limits our ability to provide an unambiguous statement on 
which portion of the parameter space in the {mg,m^o) mass plane is excluded or not, but 
nevertheless provide a very interesting estimation of the sensitivity and mass reach of the 
8 TeV analyses we reimplemented. Nevertheless, comparing with the existing literature 
enabled us to estimate the correctness of our results. First of all, comparing with Fig. 4 
of [29], we see that the shape of the monojet exclusion curve is similar and its maximum 
{mg,m^o) excluded values lie at {mg,m^o) ~ (700,700) GeV, which is quite similar to our 
own maximum which sits at {mg,m^o) ~ (750,750). In addition, the kink appearing for 
{mg,m^o) ~ (540,400) GeV is also observed in other interpretations of monojet analyses 
(see [6,29]). The underlying explanation of the kink comes from the fact that the monojet 
analysis is not, strictly speaking, based on a single hard jet. Indeed, the monojet selection 
cuts in [88] actually allow a number of jets Ajets < 3 with px > 30 GeV and |? 7 | < 2.8. Thus 
the appearance of the kink is due to the fact that, starting from this point which possesses 
a sufficient AM, additional jets fulhlling the above requirement populate the SRs, thereby 
enhancing the sensitivity. Next, comparing with the numbers given in Fig. 21 in [6], for a 
massless neutralino, gluino masses up to rrig = 1250 GeV are excluded, whereas we exclude 
up to rrig = 1300 GeV (see the broken red line in the left panel of Fig. 5). A lower mass 
limit of m^o = 550 GeV for nig = 850 GeV is also quoted to be compared with m^o ~ 600 
GeV, nig = 850 GeV on our side. Therefore our results are in fairly good agreement with 
the existing litterature. Last, there has been numerous studies, mostly in the context of 
Simplihed Dark Matter models (see e.g [107,108]), showing that the mono-jet searches can 
be sensitive to the decay width of the intermediate particles. However, since in our case the 
gluino decays through only one decay mode, its total width F is suppressed and T/nig 1, 
such that we can use safely the narrow-width approximation implicitly assumed in our SMS 
approach. 

Provided this important clarihcation, let us now discuss our results. As expected, the mono¬ 
jet analysis is the most sensitive search to probe the scenario where the gluino and neutralino 
are almost degenerate. Mass-degenerate gluino and LSP are excluded up to 750 GeV. As 
soon as AM increases, the efficiency of the monojet analysis quickly drops. This can be 
seen on the right panel of Fig. 5, where we display the two most sensitive recasted analyses 
in the small AM region. The colour coding is the same as in the left panel. The monojet 
analysis is the most effective in constraining a compressed gluino-neutralino spectrum up to 
AM ~ 50 GeV. 

For larger AM the ATLAS multijet analysis becomes the most sensitive one in constraining 
the radiative decay. At low nig the signal region dubbed as 2jl (“2 loose jets ”, more details 
about the SR dehnition can be found in the official analysis [89]) drives the exclusion curve, 
at moderate nig it is driven by SR 2jl and 3jm (“3 medium jets”) and at large nig the 
4jl (“4 loose jets”) SR dominates. Within this analysis gluino masses up to nig ~ 1300 
GeV for a massless neutralino are excluded. It appears that the CMS multijet analysis is 
less powerful than the two ATLAS ones. There are two reasons for that, one is that the 
jet binning only starts from 3 jet bins and the second is that the search was designed as a 
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discovery one, therefore its constraining power is less effective. Another comment is worth to 
be mentioned here: below the line where the tt threshold is open, interpreting the exclnsion 
contonr is highly debatable since we expect the three-body decays to take over. Therefore 
below this line onr results should be taken with caution and the official results focussing on 
the three-body decay g ^ ti are more realistic. 

In the light of these results, it would be interesting to perform a combination of at least 
the two ATLAS analyses to get a more realistic exclusion curve. However, this goes beyond 
the scope of this paper since we do not have the necessary information at our disposal to 
perform the combination. 

4 Prospects at the 13 TeV LHC 

We address the discovery and exclusion reach of the radiative gluino decay at 13 TeV in this 
section. Since the radiative decay is dominant with a small mass gap between g and Xi 2 ^ 
design a special set of cuts optimized to investigate this case. To simplify matters we select 
benchmark points with a g — Xi 2 SMS in mind, where the decay g —)■ gxi 2 is 100 %. As we 
already observed in section 2, such an assumption is well motivated from theoretical con¬ 
siderations. We choose three representative benchmark points(BP) with nearly degenerate 
masses for g and Xi ■ 

To achieve this we decouple the rest of the spectrum from this set, i.e we set the three 
generation of squarks and the sleptons to 3 TeV. Since a higgsino like scenario is favourable 
for the radiative decay, we set Mi, M 2 to 2 TeV, while /i is kept relatively light so that Xi^ X 2 
are Higgsino like. For all benchmark points, the mass gap Am = rrig—m^a is hxed to around 
10 GeV. The availability of phase space ensures that both g —)■ gxi and g —)■ gx 2 equally 
important. However, since Xi X 2 almost degenerate, the decay products of ^ are 
soft, and are most likely to go undetected. 

The relevant masses and the branching ratios (BR) are tabulated in Table 1. The benchmark 
points were generated using SUSYHIT [86]. All the benchmark points satisfy the constraints 
derived in the previous section. 



BPl 

BP2 

BP3 

M 3 

730 

820 

1100 

h 

985 

1180 

1370 

m-g 

1005 

1205 

1405 

m^o 

A1 

995 

1195 

1395 


1004 

1201 

1398 

BR{g^gxl)i%) 

99 

53 

47 

BRig^gxWo) 

1 

45 

47 


Table 1: The masses and branching ratios for the three benchmark points used for the collider 
analysis. 
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Since we focus on the nearly mass degenerate g — Xi scenario, the conventional search strategy 
is to look for a monojet + missing energy signature, where the hard ISR monojet recoils 
against the missing energy. However, at 13 TeV the number of QCD radiation jets from 
the gluino can be fairly large and carry large transverse momenta. We argue in this section, 
that a dijet + missing energy scenario (from ISR/FSR) is better suited compared to the 
traditional monojet + missing energy search strategy. We argue that the number of hard 
ISR jets can be signihcantly large at 13 TeV LHC, and including a second ISR jet helps 
considerably in improving the signal to background ratio. 






Figure 6: The distribution of the number of njets(top left), and the transverse momentum of 
the leading jet, top right) are presented in the top panel, while the transverse momentum 
of the second leading jet(p}f, bottom left) and the Mt 2 right) distribution is plotted 

in the lower panel for the signal benchmark point BP2 and the W/Z + jets backgrounds. 
The Mt 2 distribution is subject to all other cuts, as described in the text. The figures are 
normalized to unity. 

The principal backgrounds to this process are listed below, 

• QCD: The QCD dijet events are the largest source of backgrounds in terms of cross 
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section (~ 10®pb). A significant amount of missing energy is also likely to be present 
from detector mis-measurements. 

• tt + jets: The multijet + missing energy topology is fairly common for all hadronic 
top decays. It is possible that the leptons and jets are missed which can lead to a dijet 
+ missing energy topology. Since the ti + jets cross section is significantly large at 13 
TeV (~ 900 pb), even a small efficiency after cuts can result in a significantly large 
background cross section. 

• Z + jets: This constitutes the irreducible part of the background when the Z boson 
decays invisibly to a pair of neutrinos. Since the cross section is quite large ( lO^pb), 
suppressing this background to a reasonable level is of primary requirement. 

• ZZ: This background is similar to the Z + jets background, with one of the Z boson 
decaying hadronically while the other invisibly. This background is thus also one of 
the irreducible backgrounds for our case. 

• W + jets: For a leptonic decay of W, a missed lepton would mean a topology similar 
to the signal. This has to be therefore considered as a serious background for this 
topology. 

• WZ : The WZ background can be significantly large for a hadronically decaying W 
and an invisibly decaying Z. 

Parton level events for the signal and background are generated using MADGRAPH5.0.7 [96], 
and passed on to PYTHIA6 [103] for showering and hadronisation. For the signal additional 
two partons are generated at the matrix element level. For backgrounds like ti + jets, W/Z 
+ jets, additional 3 partons are generated in MADGRAPH5. The merging parameter for the 
matrix element-parton showering matching is set to 50 GeV. Since we anticipate that the 
bulk of the backgrounds are likely to originate from the tail of the pt spectrum we simulate 
W/Z + jets with a hadronic energy i77’(defined as cut of 400 GeV at the parton level in 
MADGRAPH5. The sample is then passed on to DelphesS [104] (specifically the MadAnalysisS 
tune [84] ) for detector simulation. Jets are reconstructed using Fast jet [109] with an anti- 
kr [110] algorithm and a pr threshold of 30 GeV and \r]j\ < 3. Background cross sections are 
calculated using MADGRAPH5, while the gluino pair production cross section with decoupled 
sfermions are are quoted from the official SUSY cross section working group page [111]. We 
use the MadAnalysisB [84] framework for the analysis of signal and background. 

The following cuts are applied to maximize the signal to background ratio: 

• Lepton Veto (Cl): Since the signal has no leptons, a lepton veto is imposed. Leptons 
are selected with a with transverse momentum threshold, pip > 10 GeV and |? 7 | < 2.5. 
Leptons are isolated by demanding that the scalar sum of the transverse momenta of 
all visible stable particles within a cone Ai? = 0.2, does not exceed 10 % of pip. This 
ensures that backgrounds with leptonic decays of W/Z are suppressed. 

• Jets (C2): We demand exactly two jets in the final state. From the plot of number 
of jets (njets) for the signal and background (top left Fig. 6), we observe that even 
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with a small mass gap for the signal, it is possible to obtain up to 3 jets quite easily 
from QCD radiation. The background, especially, for the W+ jets, peaks at 4 jets. 
Therefore it is convenient to design a search strategy with a dijet signature. 

• b-jet veto (C3): To suppress the large tt background, we veto b-jets. This also 
ensures that backgrounds like Zbb, Wbb are suppressed. 

• > 600,200 GeV (C4); From Fig. 6 we note that the leading ISR jet can be 
signihcantly energetic for the signal (top right. Fig. 6) and therefore we impose a hard 
cut on the transverse momenta of the leading jet. Although the second leading jet is 
not as energetic (bottom left. Fig. 6), the tail of the distribution can go up to about 
400 GeV. We therefore impose a cut of 200 GeV on the second leading jet. The hard 
cuts on the leading jets ensure that we are in a region where the missing transverse 
energy that recoil against these hard ISR jets is high. 

• Mt 2 (C5); For this study, the variable Mt 2 [112-114] is dehned as, 

MT2(ji,j2,^T) = min[maxMT(jl,flT))MT(j2,flT)] 

where the minimization is performed over all possible partitions of the missing trans¬ 
verse momentum, + and ji,j 2 , represent the two leading jets. The trans¬ 

verse mass of the system is dehned as 

M|(j, i>^) = Mj + Ml + 2(EirE^ - pip.p^). 

The value of Mt is bounded by the parent particle. It can be shown that this holds 
true for Mt 2 as well, i. e, after minimization, Mt 2 is bounded by the mass of the parent 
particle. To suppress SM backgrounds, where the only source of missing energy are the 
neutrinos, the value of the mass of the invisible particle is set to 0. SM backgrounds 
like WW, WZ, ZZ, tt, are expected to be cut off at the parent particle mass, while since 
the gluino is massive it is expected to go up to higher values. Uncorrelated/ISR jets 
however can signihcantly alter the shape of the Mt 2 distribution [115]. For low mass 
diherences, the extra ISR jets that dominate the jet conhguration are not correlated 
with the plrp, and hence the Mt 2 distribution extends well beyond the gluino mass. In 
the bottom left panel of Fig. 6 we present the Mt 2 distribution for the signal(BP2) and 
the W/Z -I- jets background subject to all the previously imposed cuts. We observe 
that the Mt 2 distribution falls sharply for the background, while the signal remains 
fairly hat. We use an optimized cut of 800 GeV to suppress the background. 

We summarize our cut hows for signal and background in Table 2. For simplicity, we only 
document the backgrounds which leave a non negligible contribution after all cuts. In Table 
2, the production cross section (G.S) for signal and background processes are tabulated in 
the third column in femtobarn. The number of simulated events (N) are noted in the 4th 
column. From the hfth column onwards the cumulative ehect of cuts are presented. 

From Table 2, we observe that the cut on the transverse momenta of the leading jets (G4), and 
Mt 2 (G5) are most ehective in suppressing W/Z-l- jets background. The other backgrounds 
(not tabulated) like QGD, ti + jets, WW,WZ,ZZ are reduced to a negligible level after the 
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mg, m^o 

G.S 

N 

G1+G2+G3 

G4 

G5 


(GeV) 

fb 





BPl 

1005,999 

314 

lOOK 

134 

4.2 

1.6 

BP2 

1205,1195 

83 

lOOK 

46 

2.3 

0.6 

BP3 

1405,1395 

24 

lOOK 

12 

0.7 

0.21 

Z -1- jets 


1.7 X 10® 

5M 

783250 

12395 

1.7 

W-b jets 


5.8 X 10® 

5M 

31921.1 

7793.1 

1.0 


Table 2: The cut flow table for the signal and background for the analysis. For simplicity 
we only provide the processes which yielda non-negligible background contribution. The cross 
section after each cut is tabulated from column 5. 


application of the Mt 2 cut. Although the signal process suffers as well as a result of the 
Mt 2 cut, the signal to background ratio is enhanced considerably. After all cuts, the total 
background cross section is reduced to 2.7 fb, while the signal cross sections vary from 1.6 
fb for benchmark PI to 0.21 fb in P3. In Table 3, we estimate the signal signihcance dehned 
as S/\fB by keeping the mass gap between the gluino and the lightest neutralino hxed at 
10 GeV for the benchmark points. 

We estimate that one can observe a signal at 5 a signihcance at 13 TeV with 30 fb“^ 
luminosity up to a gluino mass of 1 TeV while a gluino mass of 1.2 TeV can be excluded at 
the same luminosity. At 100 fb“^ luminosity the discovery reach goes up to 1.15 TeV, while 
we can exclude a gluino mass of up to 1.3 TeV. For the high luminosity run of LHC (3000 
fb“^), we can discover a gluino mass up to 1.45 TeV, excluding about 1.55 TeV. ^ 

We would like to compare our results with the standard monojet strategy at 13 TeV. However, 
since this is not available yet, a rough estimate could be obtained by a rescaling of the 
cross section from the 8 TeV result to the 13 TeV by assuming that the acceptance times 
efhciency remains the same. This has already been performed by the authors of [29] for the 
high luminosity LHC run (3000 fb“^). They predict that one can exclude a gluino mass of 
1.25 TeV with a AM of 10 GeV, although in the context where BR(^ —)■ qqx\) = 100%. On 
this evidence the dijet + missing energy strategy devised here has a clear advantage, and we 
advocate the use of this strategy for 13 TeV LHC. 

5 Conclusion 

We have performed an investigation of the MSSM region where the gluino-neutralino spec¬ 
trum is compressed and the loop-induced decay g —?■ gxfl is dominant over three-body decays. 
For small gluino-neutralino mass differences AM, this decay mode is a sensitive probe of 
compressed spectra since it does not suffer from kinematical suppression factors. Within 

^For the large mass gap scenario, we easily obtain at least 4 jets from the gluino decay and additional 
QCD radiation jets. In this case the improved CMS and ATLAS multijet searches at 13 TeV are expected 
to be quite sensitive (as we have already observed for the 8 TeV LHC in section 3). Indeed we find that our 
search strategy, which targets a dijet + missing energy signature, suffers in this case, and we only manage 
to discover a gluino mass of 1.1 TeV for a Xi mass of 1 TeV with 3000 fb“^ luminosity. 
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PI 

P2 

P3 

mg, m^o(GeV) 

1005,999 

1205,1195 

1405,1395 

^/\/R(30 fb-i) 

5.3 

2.0 

0.7 

fb-i) 

9.7 

3.7 

1.27 

^/^(3000 fb-i) 

53 

20 

7 


Table 3: The signal significances(S/y/B) for different energies and luminosities for the 
benchmark points. 


the scenario we investigated, namely with decoupled heavy scalars but with third generation 
squarks closer to the electroweak scale, the loop-induced decay of the gluino into a pure 
higgsino-like neutralino dominates over other branching fractions until the tb threshold is 
crossed, when AM increases. This occurs whatever the hierarchy between the stops and the 
sbottoms, albeit larger branching fractions are reached when 3> mj_^. On the contrary, 
the dominance of the radiative gluino decay into a bino-like neutralino or mixed bino-higgsino 
one over the other decay modes strongly depends on this hierarchy and other three-body 
decay modes start to be signihcant before the opening of the tb kinematical threshold. 
Adopting a Simplihed Model approach, where the gluino and the neutralino y? are the only 
relevant superparticles accessible at the LHC, and BR(^ —)■ gxi) = 100%, we reinterpreted 
some ATLAS and CMS 8 TeV SUSY searches to derive 95% CL exclusion contours in the 
mg — m^o plane. This decay process is particularly interesting to probe small AM. The rein¬ 
terpretation of the ATLAS and CMS analyses has been performed through the MadAnalysisS 
framework [83,84] and publicly available on its Public Analysis Database web page [95]. The 
recasted monojet analysis enables us to illustrate that, had it been used by the experimental 
collaborations to constrain the loop-induced gluino decay, this search is sensitive to degener¬ 
ate gluino-LSP scenarios. In turn, degenerate situations with gluino masses up to rUg ~ 750 
GeV could have been excluded. This surpasses the official maximal value ~ 600 GeV 
reached by the experimental collaborations using the process g —)■ qqxi- 
Then, concentrating on a pure di-jet search, we designed a dedicated search strategy to 
discover the gluino through the g —?■ gxi channel. For a mass gap hxed at AM = 10 GeV, a 
signal at 5 (T signihcance at 13 TeV with 30 fb“^ luminosity could be observed up to a gluino 
mass mg = 1 TeV, while the discovery reach goes up to 1.2 TeV. At 100 fb“^ luminosity the 
discovery reach goes up to 1.15 TeV, while we can exclude up to 1.3 TeV. A high-luminosity 
LHC could discover a gluino mass up to 1.45 TeV while excluding about 1.55 TeV. 

In addition, our Monte-Carlo investigation indicated an interesting point. Indeed, a pure di¬ 
jet search could be as sensitive as a monojet search and maybe even better. The underlying 
reason is related to the fact that a heavy gluino radiates a lot of high-p-p ISR jets, even for 
small AM. Of course this has to be taken with a grain of salt since our investigation was 
performed only a the Monte-Carlo level and only a full detector simulation could give a clear 
cut answer, but this result is nevertheless worth to be pointed out. 

The discovery of a 125 GeV Higgs boson suggests that the scalar superpartners may he in 
the multi-TeV range. Direct search of such heavy particles is impossible at the LHC. If the 
scalar mass scale fh is much less than rh < O(IO^) TeV, the gluino two-body decay is quite 
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5 Conclusion 


important in probing indirectly fh. In the end, it may well be our only probe at our disposal 
in the case of SUSY scenarios with heavy scalars and compressed SUSY, if the gluino decays 
promptly in the detector. We hope our work will motivate the experimental collaborations 
to perform a dedicated analysis to search for this decay in order to make the most of the 
LHC potential. 
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B Outflow for the ATLAS multi-jet analysis 


A Outflow for the ATLAS monojet analysis 

We give below the cut flow table for the comparison between the ATLAS monojet analysis 
and the recasted MA5 one. The monojets analysis targets the flavour-violating decay t —)■ cXi 
with a branching fraction of 100% using a monojet and c-tagged strategies. We only recasted 
the monojet strategy. The dataset corresponds to 20.3 fb“^ of integrated luminosity at 
= 8 TeV. 


Benchmark {mg,m^o) 

(200/125) GeV 

(200/195) GeV 

(250/245) GeV 


MA5 

GMS 

MA5 

GMS 

MA5 

GMS 

cut 

result 

result 

result 

result 

result 

result 

E™'®® > 80 GeV Filter 

192812.8 

181902.0 

104577.6 

103191.0 

36055.4 

48103.0 

> 100 GeV 

136257.1 

97217.0 

82619.0 

64652.0 

29096.3 

23416.0 

Event cleaning 

- 

82131.0 

- 

57566.0 

- 

21023.0 

Lepton veto 

134894.2 

81855.0 

82493.9 

57455.0 

29041.8 

20986.0 

Ajets ^ 3 

101653.7 

59315.0 

75391.5 

52491.0 

26295.2 

18985.0 

A(p{Ef^fjets) > 0.4 

95568.8 

54295.0 

70888.1 

49216.0 

24676.9 

17843.0 

PtUi) >150 GeV 

17282.8 

14220.0 

25552.0 

20910.0 

9652.1 

8183.0 

Ef^^ > 150 GeV 

10987.8 

9468.0 

21569.1 

18297.0 

8363.0 

7290.0 

Ml Signal Region 

PtUi) > 280 GeV 

2031.2 

1627.0 

4922.0 

3854.0 

2156.1 

1748.0 

Efl‘®® > 220 GeV 

1517.6 

1276.0 

4628.4 

3722.0 

2022.9 

1694.0 

M2 Signal Region 

PtUi) > 340 GeV 

858.0 

721.0 

2509.0 

1897.0 

1107.4 

882.0 

E^‘®® > 340 GeV 

344.4 

282.0 

1758.9 

1518.0 

817.5 

736.0 

M3 Signal Region 

Pr(ji) > 450 GeV 

204.3 

169.0 

773.3 

527.0 

376.1 

279.0 

E^‘®® > 450 GeV 

61.3 

64.0 

476.8 

415.0 

268.0 

230.0 


Table 4: Outflow comparison for three benchmark point of the monojet analysis targeting 
t cxj in [88], The official ATLAS numbers are taken from [116] 


B Outflow for the ATLAS multi-jet analysis 

This ATLAS analysis [89] searches for new physics in the 0 lepton -|- multi-jets (2-6 jets) 
-|- missing energy hnal state. The dataset corresponds to 20.3 fb“^ of integrated 

luminosity at = 8 TeV. In the context of supersymmetry, the analysis targets gluino pair 
production {gg), squark pair production (gg), and squark gluino production(gg). For each 
of the production modes the following cases are investigated; 
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B Outflow for the ATLAS multi-jet analysis 


decay to LSP 

2jm 

(qq) 

2jt 

(qq) 

3j 

iM) 

4jl- 

{m) 

4jl 

{E) 

4jt 

{9~9) 


direct 

direct 

direct 

direct 

direct 

direct 

f^glg 

475 

1000 

1612 

400 

800 

1425 

^xt 

- 

- 

- 

- 

- 

- 

m^o 

_dj_ 

425 

100 

337 

250 

650 

75 


Table 5: The benchmark scenarios used for the valiaion of the 10 signal regions. 

For the direct sguark pair or gluino pair production, the rest of the spectrum is de¬ 
coupled. For sguark-gluino production, the gluino and sguark masses are assumed 
to he degenerate. All masses are in GeV. The SLFIA files were obtained from 

http: //hepdata. cedar, ac. uk/view/ins1298722 . 



2jm 

2jt 

3j 

4jl- 

4jl 

4jt 


15 

15 

- 

10 

10 

- 


- 

- 

0.3 

- 

- 

- 

Meff(incl) 

1200 

1600 

2200 

700 

1000 

2200 


Table 6: The signal region specific cuts. All energy units are in GeV. 


• Gluino pair production : 

1. Direct: g —?■ ggy?. 

2. One step: g —)■ qq'xf followed by xf ^ Wxi 

• Squark pair production : 

1. Direct : q -7 qxi- 

2. One step : q —?■ q'x^ followed by xf ^ Wxi 

• Squark gluino production : 

1. Direct : g ggy?, q qXi- 

For all of the above cases, all branching ratios are assumed to be 100 % and the rest of the 
spectrum is decoupled from this set. 

In the following tables we summarize the benchmark points used (Table 5), the cuts used 
for different signal regions(up to 4 jets, Table 6) , and hnally the comparison between the 
ATLAS multi-jet analysis and the recasted MA5 one (Tables 7-8). Further details about the 
validation can be found on the web page [95]. 
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B Outflow for the ATLAS multi-jet analysis 


Benchmark/SR 

2jm 


2jt(gg) 

3j(^g) 


MA5 

ATLAS 

MA5 

ATLAS 

MA5 

ATLAS 

cut 

result 

result 

result 

result 

result 

result 

> 160, 

PT{ji,j2) > 130,60 

1656.1 

1781.2 

62.1 

61.6 

18.8 

18.6 

PtUs) > 60 

- 

- 

- 

- 

15.1 

14.8 

A0(i,, E”-) > 0.4 

1295.9 

1462.7 

56.9 

55.7 

13.3 

12.9 


449.1 

566.1 

40.1 

38.5 

- 

- 


- 

- 

- 

~ 

10.1 

9.6 

Meff(incl) 

122.2 

102.4 

23.8 

21.7 

6.2 

5.9 


Table 7: Outflows for signal regions 2jm,2jt and 3j, compared to the official ATLAS results 
documented in [89], All energy units are in GeV. 


Benchmark/SR 

cut 

4jl-( 

ATLAS 

result 

qq) 

MA5 

result 

4jl(5 

ATLAS 

result 

9 ) 

MA5 

result 

4jt(c 

ATLAS 

result 

Eg) 

MA5 

result 

Etf^^ >160, 

PtO'i, 72 ) > 130,60 

16135.8 

15097 

634.6 

679.0 

13.2 

12.7 

PtUsOa) > 60 

2331 

2112 

211.4 

185.7 

12.0 

12.0 

> 0.4 

A0(j4, > 0.2 

1813.7 

1723.0 

154.6 

144.9 

8.4 

8.9 


1009 

943 

98.7 

84.4 

- 

- 


- 

- 

- 

- 

4.8 

5.5 

Meff(incl) 

884 

843 

39.5 

41.5 

2.5 

2.9 


Table 8: Outflows for signal regions 4jl-,4jl and fjt, compared to the official ATLAS results 
documented in [89], 
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